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ABSTRACT 

The apparent age and mass of a stellar cluster can be strongly affected by stochastic sampling of the 
stellar initial mass function, when inferred from the integrated color of low mass clusters (< 10 4 M ). 
We use simulated star clusters to show that these effects are minimized when the brightest, rapidly 
evolving stars in a cluster can be resolved, and the light of the fainter, more numerous unresolved 
stars can be analyzed separately. When comparing the light from the less luminous cluster members 
to models of unresolved light, more accurate age estimates can be obtained than when analyzing the 
integrated light from the entire cluster under the assumption that the initial mass function is fully 
populated. We show the success of this technique first using simulated clusters, and then with a stellar 
cluster in M31. This method represents one way of accounting for the discrete, stochastic sampling 
of the stellar initial mass function in less massive clusters and can be leveraged in studies of clusters 
throughout the Local Group and other nearby galaxies. 

Subject headings: galaxies: stellar content 



1. INTRODUCTION 

Stellar clusters are used as tools in the study of star 
formation, stellar evolution, and galactic evolution be- 
cause they contain numerous stars of the same age, dis- 
tance, and metallicity. Their scientific utility rests on 
our ability to accurately measure their ages, masses, and 
metallicities. 

Resolved observations of individual stars provide the 
most direct way of measuring the properties of clusters. 
By resolving individual stars, one can analyze the re- 
sulting color-magnitude diagrams (CMDs) by comparing 
them to isochrones, allowing one to derive the age, ex- 
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20071. An additional advantage is that one can obtain 
an estimate of extinction that is independent of model 
colors. However, this approach is limited by the angular 
resolution of the available images. If a cluster is too dis- 
tant, there will not be enough individual stars resolved 
to use isochrone fitting. The distance limit for analyz- 
ing clusters through resolved stars varies with the surface 
brightness of the cluster and the resolution of the tele- 
scope, but in practice isochrone fitting has rarely been 
used beyond a few Mpc, even with the Hubble Space 
Telescope (HST). 

If the stars in a cluster are unresolved, then the prop- 
erties of the cluster can only be de rived from th e inte- 
grated light of the entire cluster (e.g., Larsen 2009). Such 
studies have less stringent requirements tor angular reso- 
lution, and can be carried out at larger distances (up to 
~ 50 Mpc; Adamo et al.||2010" ). Thus, analysis of the in- 
tegrated light of unresolved star clusters will remain the 
method of choice for deriving ages and masses of clusters 
in more distant galaxies, which span the largest possible 
range of environments and star formation modes. 

There is growing awareness of the potential limitations 
of unresolved cluster studies. One of the most severe is 
stochastic sampling of the stellar initial mass function 
(IMF), which strongly affects the integrated colors and 
fluxes of low mass clusters. Discrete sampling of the IMF 
affects the present day mass function, altering a cluster's 
color and luminosity from what would be expected if the 
IMF were continuously populated up to the highest pos- 
sible stellar mass, and causes potentially rapid changes in 
color and magnitude throughout the lives of the cluster 
( Piskunov et al.||2011[ ). Clusters with masses less than a 
few 10 4 Mq are subject to the effects of small number 
statistics at the upper end of their mass function, and 
thus their light can be dominated by just a few post-main 
sequence stars. A small variation in their number or evo- 
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lutionary state can cause two otherwise similar clusters 
to have drastically different integrated colors (up to five 
magnitude s ) (e.g., Barbaro fc Bertelli||1977[ |Girardi fc| 
Bical[l993l |Girardi et alJ|l(WB| pantos k Frogel||l9^7[ 
Brocato et al.||1999| |Langon fc Mouhcine||2000l |Bruzual| 
2002| |Anders et al |I2004| |Raimondo et al |2005] |Ueryino 
& Luridiana 2006; Dcveikis et al. 2008; Popescu & Han- 
son|2010b[ see also Figure [l] in ij2|). This spread in colors 



leads to large systematic errors when deriving the clus- 
ters' properties. Stochastic effects have been shown to 
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iNew techniques nave been developed, to more accu- 
rately age-date low mass clusters, by taking stochastic 
sampling into account. These improvements are impera- 
tive, given that current surveys are delving more deeply 
into the low mass regime. Sp ec ifically, recent studies 



such as De veikis et al. rt2008L Fouesnea u & Langon 



(2010), and Popescu fe Hanson| pl)10b|a ) have moved 



away from the traditional population synthesis mod 
els towards discrete population models, where the mass 
function is explicitly considered as a discrete distribution 
of stars. The main ingredients of such techniques are 
large collections of discretely sampl ed synthetic cluster s 
from Monte-Carlo simulations (e.g., |da Silva et al.|2012 ), 
which can be compared with the observed photometry. 
Alternatively, individual filters ca n be down- weighted t o 
reflect the degree of stochasticity ( Mafz Apellaniz[ 2009 ) . 

The growing body of literature on stochastic iiuctu- 
ations of clusters has clearly shown that these fluctua- 
tions are due to the random presence of massive upper 
main sequence and evolving stars within a population. 
There is also evidence that the flux from the main se- 
quence is f airly stable in compa rison (within a few tens 



of percent; Langon et al. 



20081. As the quality of the 



observations has improved, we are currently able to re 
solve the most massive cluster stars out to distances of 
about 4-5 Mpc. In many cases, one can explicitly iden- 
tify the tiny subpopulation of rapidly evolving luminous 
stars, and exclude this small stochastically sampled pop- 
ulation from analysis of the integrated light. Removing 
the post-main sequence stars from the analysis reduces 
the effect of stochasticity, offering a promising way for- 
ward for deriving accurate ages of low mass clusters. 

In this paper, we explore a method for analyzing only 
the unresolved part of a young cluster's light, which in 
the ideal case includes only stars up through the main 
sequence turnoff. This approach minimizes some of the 
troublesome effects related to the stochastically sampled 
upper end of a cluster's stellar mass function. Since 
younger clusters are more susceptible to stochastic fluc- 
tuations and are found in abundance in nearby galaxies, 
we focus on minimizing errors for the younger clusters. 
We demonstrate the effectiveness of this method using 
synthetic clusters. We then apply this method to a clus- 
ter in M31 from the Panchromatic Hub ble Andromeda 
Treasury (PHAT; [Dalcanton et alT][2012| . 

This paper is organized as follows: In ^2] we demon- 
strate the problem of stochastic fluctuations in the inte- 
grated colors of clusters. In Sj3]we describe our method 
of using only unresolved light to study partially resolved 



clusters. We show the results of a simultaneous deriva- 
tion of age, mass, and extinction for simulated clusters 
in <HJ In ^5] we test our method on a cluster in the 
PHAT survey which has independent measurements of 
its parameters. We discuss practical aspects of using 
this method in ^6] and state our conclusions in £j7J 

2. STOCHASTIC FLUCTUATIONS 

A cluster's light is dominated by the brightest, most 
massive stars. A small stochastic difference in the num- 
ber or evolutionary state of these stars can cause two 
otherwise similar clusters to have drastically different in- 
tegrated colors. This effect is very pronounced in young 
clusters (< 100 Myr), whose flux may be dominated by 
a small number of bright red giants or supergiants that 
bias the integrated color towards the red. To illustrate 
the effects of stochasticity in a young cluster, we created 
a synthetic example that shows the dramatic impact that 
one evolved star may have on a c luster's integrated light , 
similar to the examples shown in Santos & Frogel ( 1997). 
For consistency with the PHAT survey, we will consider 
filters used in this survey, in particular F336W, F475W, 
F814W, and F160W. 
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Simulated color-magnitude diagrams for two 10 Myr 
old, 10 a Mq clusters at solar metallicity. Due to stochastic sam- 
pling of the initial mass function, the cluster in the left panel does 
not contain any evolved stars, while the cluster on the right has 
one evolved red giant. Integrated colors are included for each panel, 
along with a 10 Myr isochrone, and a horizontal line showing the 
magnitude of the main sequence turnoff. 

Figure [l] shows two realizations of a 10 Myr old clus- 
ter with a mass of 10 3 M F)l generated using the pro- 
gra m F ake (Dolphin 2002J, which is described in detail 
in §3.1[ Both of these clusters are generated assuming 
the same underlying isochrone, total mass, and initial 
mass function. The thickness of the main sequence is 
due to binarity , assuming a binary fraction of 0.35 and a 
Salpeter| ( |1955[ ) slope for the binary mass distribution. In 
the realization in the left panel, the sampled population 
does not include any rapidly evolving post-main sequence 
stars. The integrated (F336W-F814W) color is -1.56, 
and is dominated by the massive blue stars at the top 
of the main sequence. In the right panel, however, the 
cluster (again, with the same mass and age) happens to 
have one very massive, bright red giant, due to stochastic 
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Fig. 2. — Examples of the va riation in cluster resolution with distance. The left panel shows the M31 cluster PC1017, at a distance of 
0.785 Mpc jjohnson et al.|2012| . The middle and right panels show image simulations of this cluster at roughly twice and four times the 
distance of M3T (1.57 and 3.14 M 
distances. 



Mpc). The image is in the F814W filter, and has been smoothed and rebinned to simulate the two further 



sampling of the initial mass function. This bright evolved 
star now dominates the total integrated color of this clus- 
ter, which is now +2.77 - a difference of more than four 
magnitudes. A 10 Myr old cluster could therefore po- 
tentially have the same integrated color as a 10 Gyr old 
cluster due to one evolved giant star. Although the effect 
is smaller, one might also expect stochastic variations in 
the blue flux due to variations in the sparsely populated 
upper main sequence. This can also be seen in Figure [l] 
where there is a 1 magnitude difference in the luminos- 
ity of the brightest main sequence star between the two 
cluster realizations. These effects make age and mass 
determinations based on integrated color extremely un- 
certain. 

While Figure [l] shows large differences in the two clus- 
ter realizations brightward of Mpsuw = —2, faintward 
of this limit the two clusters are populated quite simi- 
larly. This similarity suggests that if the single red giant 
could be excluded, the color and luminosity of the rest 
of the stellar population would be better behaved, and 
would share the same blue colors expected for a young 
cluster. This could in practice be done if the brightest 
one or several stars that correspond to the red giant and 
supergiant phases are resolved. 

With HST's high resolution imaging, there are a num- 
ber of large nearby galaxies whose clusters are partially 
resolved into stars. Figure [2] demonstrates the variation 
in cluster resolution with distance. The left panel shows 
the M31 cluste r PC1017, at a distance of 0.78 Mpc ( [John- 
son et al.|2012|) . The middle and right panels show image 



simulations ol this cluster at roughly twice and four times 
the distance of M31 (1.57 and 3.14 Mpc). Even at these 
distances, the cluster's brightest stars can be resolved 
(in particular, the red evolved stars), while the majority 
of the more numerous main sequence stars remain un- 
resolved. These fainter stars are blended together and 
fall below the detection limit, particularly in the center 
of the cluster, where the crowding and background level 
are high. Many older clusters at larger distances may 
also fall in this regime, because evaporation and cluster 
dissolution lead them to have lower surface densities, al- 
lowin g individual bright stars to be detected ( |Gieles et al 
20111. 



3. METHOD 

3.1. Simulated Clusters 

To demonstrate the potential of analyzing the unre- 
solved flux in a partially resolved cluster (rather than 



its fully integrated counterpart), we start by building a 
set of unresolved light models. This includes choosing a 
val ue fo r the magnitude cutoff, My lm , which is described 
in {3.2 Then, based on this set of models, we show our 
method using simulated clusters and an M31 test cluster 
from the PHAT data set. 

We utilized a variety of simulated clusters to test the 
degree to which stochastic effects arc suppressed by re- 
stricting the light to the unresolved portion. We gen- 
erated CMDs for simulated clusters using the program 
Fak e, which is par t of the CMD analysis suite MATCH 
(see ~ 
synt 



Dolphin 
retic ste. 



2002 1 . This program generates CMDs for 



ar populations usi ng theoretical isochro nes 



We used the Padova models (Marigo et al. 2008) for 



the isochrone set, with updated AGB tracks from Gi 



rardi et al. 



20101 along with a Salpeter initial mass 
peter 1955). The Salpeter IMF produces 



function 

more low-mass sta rs and a higher mass-to-light ratio than 
the Kroupa IMF (|Kroupa||2001|). The simulations were 
done using HST filters, since we plan to apply this tech- 
nique to M31 clusters observed with HST. We chose to 
study clusters at four fiducial masses: M c i / M Q = 10 3 , 
10 4 , 10 5 , and 10 6 . This range is representative of the 
various masses of clusters we expect to be able to de- 
tect in nearby galaxies with HST. At the lower mass 
end, the clusters experience a high level of stochastic- 
ity, while at the high mass end, clusters fully populate 
their isochrones, and suffer few color and luminosity vari- 
ations due to stochasticity. The simulated clusters' ages 
t ranged from log(i/yr) = 7 to log(i/yr) = 10, in incre- 
ments of A log(i/yr) = 0.3 — 0.4. Thirty clusters at each 
age and mass were produced at solar metallicity, and an- 
other thirty at subsolar metallicity (Z = 0.1 Z Q ). Solar 
metallicity clusters are mainly used throughout this pa- 
per since this is appropriate for young clusters in nearby 
large spirals like M31. 

Fake generates clusters by randomly drawing masses 
from the initial mass function until the desired cluster 
mass is reached. However, if a massive star is drawn 
which will put the cluster over the given limit, this star 
is discarded in favor of lower mass stars. This approach 
leads to a bias in the resulting cluster's mass function, 
which will include more low mass main sequence stars. 
This is a known effect of the cluster simulation, and 
should be kept in mind when analyzing resulting clus- 
ter fluxes and colors. However, we are not studying the 
clusters' mass functions, and this does not affect the self- 
consistent analysis on synthetic models that we conduct 
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Fig. 3. — Color-magnitude diagrams in two of the PHAT survey filters for solar metallicity clusters at four fiducial ages. Horizontal 
lines indicate three possible values for -Mi; m : —2.0, —3.0, and —4.0. Stars fainter than these lines are considered part of the unresolved 
component. The black clusters have no extinction, while the red clusters have Ay = 2.0 mag. 



in this paper. Alternate methods exist to sampl e the 
cluster mass function, including those found in flMaiz 
Apellaniz|l2009l JWeidner fc Kroupa|2006l [20041 [Popescu 
fc Hanson|2009 ) , some of which produce better sampling 
of the mass function of very low mass clusters. However, 
the differences between our method and a proper filling 
of the cluster mass function are small for the mass range 
we are considering (M ci > 10 3 M ). 

3.2. Selection of Magnitude Cutoff 

To study the broad behavior of partially resolved clus- 
ters, we chose a simple magnitude cutoff Mn m that de- 
fines the threshold below which the cluster stars are 
treated as unresolved. This threshold does not neces- 
sarily have to be assigned to the actual magnitude limit 
of the data, and should be far from the completeness 
limit. Such an approach can compensate for the strong 
radial variation in the detection limit, which tends to 
be significantly brighter in a cluster's crowded inner re- 
gions. Ideally, we would like our value of M um to be at 
a level that excludes all the luminous evolved stars, and 
that leaves stars up through the main sequence turnoff. 
If our Mn m is too faint, the main sequence turnoff will 
not be included, thus excluding the source of the age in- 
formation. If our imposed My lm is too bright, we will 
include more post-main sequence stars, which will not 
be as effective in reducing stochasticity. In this section 
we examine the best choice of Mu m for the PHAT data 
set. A similar exercise should be done when applying 



this method to other data sets. 

Figure [3] shows the CMDs of simulated clusters at four 
fiducial ages, along with three values of Mn m : —2.0, —3.0, 
and —4.0. These values of M\i m were chosen to be above 
the detection thresholds of individual stars in M31. The 
black points are for stars in a cluster with no extinction, 
and the r ed points are for A v = 2.0 mag (assuming R v — 



3.1 and a Cardelli et al. (1989| extinction law, which we 
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use throughout this paper). The number of stars being 
counted as part of the unresolved light clearly depends 
upon the age, extinction, and value of My lm . Clusters 
younger than 25 Myr have all evolved stars cut off, re- 
gardless of Mii m and extinction. For heavily extincted 
50 Myr old clusters, some evolved stars remain below 
the less stringent values of M\i m , and are therefore not 
removed from the unresolved component. For clusters 
with little to no extinction, the age at which evolved 
stars begin to fall below the brightest M\i m is 100 Myr. 
We therefore expect that this method will offer the most 
improvements when estimating the ages of younger clus- 
ters, which have F814W turnoff magnitudes of —2 and 



3 for ages of 35 Myr and 18 Myr, respectively ( Girardi 
et al.|2010| . Since younger clusters also suffer more from 



stochastic effects, we therefore focus on minimizing er- 
rors for younger clusters. The tests described in jQshow 
reduced errors for My lm — — 2 and —3. For our analysis, 
since Mu m = —2 and —3 are both effective at reduc- 
ing stochastic effects in the young clusters, we choose to 
use the brighter My lm = — 3 to ensure that all the stars 
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Fig. 4. — Absolute magnitudes in the F814W filter as a function of age for simulated solar metallicity clusters at four fiducial masses. 
The left panel shows the total integrated flux, while the right panel shows only the unresolved flux contained in stars fainter than the 
magnitude cutoff at Mij m = —3. 30 clusters are plotted for each age. The clusters in the right panel exhibit much less spread in magnitude 
at a given age and mass, since the most luminous, rapidly evolving stars are excluded. 

again restricting the unresolved light to stars fainter than 
Mim = -3 in the F814W filter. 

Figure [5] shows that the spread in colors becomes small 
(less than 0.5 magnitudes) as the cluster mass approaches 
fO 6 M and the cluster's mass function becomes less 
subject to stochastic sampling. Even for clusters of 10 5 
M , the colors of the youngest clusters span a range of 
0.9 mag, which can lead to significant errors (more than 
100 Myr) when deriving ages from these colors. These 
uncertainties can cause significant problems in studies of 
more distant galaxies, for which clusters like these are 
the only clusters detected, and for which no individual 
stars can be resolved. 

At even lower masses (< 10 4 M@), stochasticity has 
an even larger effect, for the younger clusters especially. 
The youngest clusters exhibit a range in colors from 
(F336W-F814W) = -2 to +3, a spread of five magni- 
tudes. At red colors, the 10 3 M & clusters have only one 
or two red giants or supergiants that are responsible for 
biasing the overall color. The bluer clusters have no red 
(super)giants, such that their light is dominated by the 
and B stars on the main sequence. This again causes 
the bimodal distribution in color seen for the young low 
mass clusters. 

In contrast, when only the unresolved light is consid- 
ered (right panel of Figure [5]) , the variation in color de- 
creases to less than one magnitude for clusters with ages 
less than 100 Myr, even for the least massive clusters. 
Thus, using only unresolved light reduces variation in 
color by four magnitudes. For older clusters, the vari- 
ation in color for the least massive clusters increases to 
almost two magnitudes, since not as many of the red 
evolved stars' light are excluded. However, in practice, 
old low mass clusters are rare since most have already 
dispersed. Ignoring the older (> 100 Myr), least massive 
clusters, the color variations are on the order of one mag- 
nitude or less for all clusters of the same age and mass 
when only their unresolved light is considered. 

In addition to much smaller dispersion, the unresolved 
colors show several differences when compared to the tra- 
ditional simple stellar population (SSP) models (shown 



brighter than this will be well resolved. 

The clusters in the PHAT survey are an excellent ap- 
plication for this method since they are observed in six 
filters with HST, and the completeness is high even at 
M 8 i4 = 0. At larger distances, this method can be used 
in galaxies that meet this resolution limit, but do no have 
enough resolved cluster stars to do isochronc fitting. It is 
also important to note that the cluster photometry needs 
to be complete at or brighter than Mn m . This condition 
may be difficult to achieve for very dense clusters. 

3.3. Simulating the Unresolved Light 

Figure [4] shows the relationship between a cluster's ab- 
solute magnitude and mass, for the integrated F814W 
magnitude of both the entire cluster (left panel) and the 
unresolved light below the magnitude threshold of Mu m 
= —3 (right panel). For a given cluster mass, the vari- 
ance in absolute magnitude depends on the cluster's age 
and number of evolved stars. For the youngest, least 
massive clusters, the spread is almost five magnitudes, 
while the unresolved light spans less than 1.2 magni- 
tudes. This large magnitude dispersion is consi stent with 
previous work on stochastic fluctuations (e.g., Santos fc 
Frogelj [19971 |Piskunov et~aTl [20091 |Popescu & H" anson 



2010bp . If a young low mass cluster contains only one 
or two bright, evolved giants or supergiants, that clus- 
ter's brightness can increase by several magnitudes, and 
it may appear as a more massive cluster. In contrast, the 
clusters with no evolved stars have much fainter magni- 
tudes. This causes the bimodal distribution seen for the 
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The stochastic effects visible in Figures [T] and H] make 
age determinations using integrated colors challenging. 
Theoretically, it is possible to obtain much better results 
by looking only at colors derived from the unresolved 
light. To test this, in Figure [5] we compare the colors de- 
rived from either the integrated light (left) or the unre- 
solved light (right) to the simulated clusters' input ages, 
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Fig. 5. — Colors as a function of age for simulated solar metallicity clusters at four fiducial masses. The left panel shows the total 
integrated color, while the right panel shows only the unresolved color (excluding stars with F814W magnitudes brighter than Mr lm = —3). 
30 clusters are plotted for each age. The break in the right panel at log(t) > 8.4 is due to Mi; m being above the main sequence turnoff for 
older clusters, such that some evolved stars are included in the unresolved component. For clarity, all the 10 6 Mq clusters were shifted to 
the left by 0.06 dex, the 10 s Mq clusters were shifted to the left by 0.02 dex, the 10 4 Mq clusters were shifted to the right by 0.02 dex, 
and the 10 3 Mq clu sters were shifted to the right by 0.06 d ex. The black curve shows the colors predicted from continuously populated 
Padova SSP models (Marigo et al.|2008| |Girardi et al.|2010} . 



as the black curve) . The unresolved colors are in general 
bluer than the colors given by the SSP models, since not 
all of the evolved, red (super)giant stars' light is included. 
The unresolved colors also show a smooth variation for 
ages less than log(i/yr) = 8, while there is a modest dis- 
continuity in unresolved colors seen between log(t/yr) = 
8 and 8.4. This effect is a result of where Mum is in rela- 
tion to the main sequence turnoff (see Figure |3| and the 
age step size used. For clusters younger than log(i/yr) = 
8.4, the unresolved component excludes all evolved stars, 
while for older clusters, the choice of Mi; m leaves some 
evolved stars in the unresolved component. The colors 
of the older clusters, therefore, match more closely with 
those predicted by the SSP models. 

Similar tests conducted at a metallicity of Z = 0.1 Z Q 
yielded qualitatively similar results to the solar metallic- 
ity case presented above. 

4. TESTS ON SIMULATED CLUSTERS 

In practice, determining a cluster's properties is often 
done through a simultaneous derivation of age, mass, and 
extinction determined using y 2 m inimization compared 



to a set of fiducial m odels, (e.g., Pasquali et al. 2003 



Hancock et al.||2008[ ). The value of x\ f° r comparing 



data to a set of fe model parameters is given by 
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where Di is the data flux in the ith filter, Oik is the 
fcth model flux in the ith filter, Oi is the variance in the 
ith filter, M is the mass, and the summation is over the 
available filters. The parameters of the model k are then 
varied until x\ is at a global minimum. 

A x 2 minimization was used to simultaneously recover 
the age and extinction for a set of test clusters to a set 
of model clusters, while optimizing for the mass, using 



only their unresolved flux. The test cluste rs w ere drawn 
from the clusters used in the analysis in §3.3[ with ad- 
ditional clusters between the ages 10 My r and 100 Myr 
added to better probe the rapidly evolving stellar pop- 
ulations of young clusters. Thirty clusters at each age 
were chosen as test clusters, for a total of 540 test clus- 
ters. The test clusters' masses ranged from M c ; / Mq 
= 10 3 - 10 6 . We looked at four HST filters that are be- 
ing used in the PHAT survey: F336W, F475W, F814W, 
and F160W. Each test cluster was then given a random 
value of foreground extinction between Ay — to Ay — 
3.0 mag; no differential extinction was included for this 
initial test. A small amount of random Gaussian noise 
(<r = 0.05 mag) was also added to these clusters to simu- 
late observational errors, which we assume are Gaussian. 
This noise value was used for the variance, <Ji. through 
the conversion <jfi ux = flux x (1 — l[p - 4o ' mti 9). 

The set of clusters that were used as models were taken 
from the most massive clusters at each age. These model 
clusters were massive enough for their main sequence to 
be considered fully populated (> 10 5 Mq). Altogether, 
we used model clusters at eighteen different ages from 
log(t/yr) = 7 to 10. 

The models used were normalized as the unresolved 
flux per unit mass and, assuming that light scales as 
mass, the test clusters' unresolved fluxes were scaled to 
match the model clusters' fluxes. Therefore, mass is a 
parameter solved for during the minimization by calcu- 
lating the scale factor between the fluxes of the test clus- 
ters and the model clusters. Extinction was recovered 
by including these same models at 50 values of extinc- 
tion between Ay — to 3.5 mag. 

Figure [6] shows the recovered ages for four different val- 
ues of A/iim, color coded by the difference in recovered 
and input extinction. For Mu m = —2, the bias in re- 
covered ages is close to 0, and the average error is 0.09 
dex. The average error in mass is 0.05 dex, and there 
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Fig. 6. — Recovered vs input age of the simulated test clusters, color coded by the difference in recovered vs input extinction, for four 
values of M\i m in the F814W filter. To clearly see the distribution of recovered values, Gaussian noise with a dispersion of a = 0.05 dex in 
each direction was added in the upper panels. The bottom panels show the median for the residuals in log(t) for all clusters at each input 
age, along with the 16 th and 84 th percentile of the residuals. 



is a slight bias toward underestimating the mass. The 
average error in Ay is 0.12 mag. The errors were cal- 
culated by taking the mean of the absolute differences 
in input — output. These values represent a lower limit, 
since the calculated errors are in many cases dominated 
by the grid size. 54% of the clusters' ages are recovered 
to better than the resolution of the grid of models, and 
96% are recovered within 0.5 dex. Mu m = —3 has simi- 
lar errors, while Mu m = — 1 has a slightly larger error in 
age. M um = —4 does not offer as large an improvement, 
and shows a large number of outliers. With this brighter 
limit, evolved stars are still included in the unresolved 
light at most ages, so the stochastic effects are not as 
reduced. For fainter magnitude limits (Mi; m > — 1), the 
age sensitivity is reduced since all of the evolved stars and 
some of the bright main sequence stars are being removed 
for clusters younger than 50 Myr with low extinction, as 
seen in Figure |3[ For the least massive clusters (<10 3 
M ), the average error in age is 0.17 dex. The feature 
seen at a recovered age of log(t/yr) = 8 consists mostly 
of low mass clusters that are not well recovered. This 
occurs because, due to the shape of the continuous SSP 
models, when the minimization routine searches for the 
closest match to the models, certain age models are more 
attractive to a wider variety of clusters, which was was 
pointed out in Section 4.3 of |Fouesneau fc Langon|[2010} 
Figure [7] compares the results of the x z minimization 
for Mii m = — 3 (top panel), with the results for tradi- 
tional fitting to integrated light models (bottom panel). 
The left plots show recovered ages and the right plots 



are the recovered masses, color coded by the difference 
in recovered and input extinction. In addition to the fea- 
ture seen at a recovered age of log(t/yr) = 8, explained 
above, there is an additional feature at log(i/yr) = 10. 
This feature results from the finite range of model proper- 
ties. For clusters that are redder than the age-extinction 
range covered by the models, they are always fit to the 
oldest age in the models. This effect can also be seen 



in |Fouesneau et al.]|2012 The large feature at a recov- 
ered age of log(i/yr) = 10 is not seen in the results from 
using the unresolved light (top panel). The feature at 
log(i/yr) = 8 is also significantly reduced compared to 
the integrated light results. Due to the monotonic varia- 
tion in unresolved color, as seen in Figure [5l many of the 
sources of color degeneracies in integratea light are not 
present for the unresolved light. This leads to greater 
accuracy in deriving the cluster ages and masses. 

For the fits to integrated light, there is a bias of 0.09 
dex toward older ages, and the average error for all test 
clusters is 0.26 dex in age. 36% of the clusters' ages are 
recovered to better than the resolution of the grid of mod- 
els, and 83% are recovered within 0.5 dex. Some clusters' 
recovered ages can be as much as 3 dex away from their 
input age. The mass estimates are more robust and have 
an average error of 0.19 dex. However, significant scat- 
ter exists for the integrated light fits for the lower mass 
clusters, which violate the assumption that mass scales 
monotonically with luminosity (see Figure El), as we have 
implicitly assumed in our fitting process. The average er- 
ror for Ay is 0.45 mag, and recovered values of Ay are 
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Fig. 7. — Recovered vs input age (left plots) and mass (right plots) of the simulated test clusters, color coded by the variation in recovered 



extinction. The top row shows the results from using unresolved light with Mn n 



-3, while the bottom row shows the results from using 



integrated light. We added Gaussian noise with a dispersion of a = 0.05 dex in each direction for clarity. The bottom panels show the 
median (dot) along with the 16 th and 84 th percentile (bar) values of the residuals in log(t) and log(Af), respectively. At each age time 
step, the age residuals are shown for each input mass, where the more massive clusters' residuals are shifted slightly to the left and the less 
massive clusters are shifted slightly to the right. The mass residuals are shown for four age groups (limits of 7.0, 7.5, 8.0, 8.5, 9.0), where 
the older clusters' residuals are shifted slightly to the left and the younger clusters' residuals are shifted slightly to the right. 



biased by 0.16 mag toward lower extinctions. This ex- 
tinction bias explains the age bias, since young, highly 
extincted clusters are being recovered as older, less ex- 
tincted clusters. When considering only clusters whose 
masses are < 10 3 M Q , the average error in age is 0.51 dex, 
considerably higher than the error for low mass clusters 
when using unresolved light. 

In this analysis we explored clusters with a range of 
extinction up to Ay = 3.0 mag, which is appropriate for 
most M31 clusters. If independent constraints on Ay 



are available, the age-extinction degeneracy is lifted, and 
the errors will be reduced even further. For example, for 
-Mlim = —2, if Ay is known to within 1 magnitude, the 
error in age decreases from 0.09 to 0.03 dex. Also, us- 
ing additional filters would i mprove the accuracy of the 
age and mass determination (Mafz Apellaniz 2009 ), while 
providing greater constraints on extinction. The metal- 
licity of all of our clusters was fixed at solar, however 
this analysis could be extended to varying metallicities. 
Allowing the metallicity to vary would cause further de- 
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(F475W - F814W) 



Fig. 8. — Test cluster PC1017 in M31, as observed with the PHAT survey jjohnson et al.|2012| . Left panel is the image shown with the 
photometric aperture. Right panel shows the color-magnitude diagram for this same cluster. Trie best-fit isochrone (20 Myr, Ay = 0.9) 
from fitting the resolved stars is shown in black. The stars excluded in the unresolved fitting are shown in red. The unresolved method 
gives an age of 40 Myr and Ay = 0.42. The blue isochrone (80 Myr, Ay = 0.28 mag) that was obtained using traditional integrated light 
fitting is a poor fit to the resolved stars. 

generacies with age and extinction. A full characteriza- 
tion of errors and systematics is beyond the scope of this 
paper and is postponed to a further publication. 

5. APPLICATION TO A REAL CLUSTER 

In this section we extend the % 2 fitting to an observed 
cluste r from the PHAT cluster catalog (Johnson et al. 
20121) PC1017. Figure [|j shows the CMD for cluster 
FC1017 which we use as a test cluster to demonstrate 
our unresolved method. 

The challenge of using a real cluster to test this method 
is that we do not know its intrinsic properties. We chose 
this cluster as our test cluster because there are good 
indepe ndent estimates of its parameters. Caldwell et al. 
(2009) used spectroscopy to age-date this cluster as 33 
Myr, with a factor of 2 uncertainty. Its mass was found 
to be log(M ci / M Q ) = 3.88, and A v was 0.68 mag. 
Stochasticity can however affect these spectroscopic esti- 
mates as wel l. Another prediction given by the discrete 
models from |Fouesneau fc L ancpn 2010 gives an age of 
41 Myr, log(M ci / M Q ) = 3.9, and Ay = 0.3 mag. Ad- 
ditionally, since this cluster has many resolved stars, we 
were able to make determinations by isochrone matching 
the best fit age and extinctio n with the CMD analysis 
suite MATCH ( |Dolphin|[2002| . This gives an age of 20 
Myr (log(Vyr) = 7.3 +7- 6.1), log(M ci ) = 3.8 +/- 0.3, 
and an Av = 0.9 +/- 0.15 mag. 

We measure the unresolved flux for the cluster using 
a three step process. Due to large photometric errors 
in F160W, we use only F336W, F475W, and F814W. 
For each of the six images, we sum the light within the 
photometric aperture (2.19") and within a background 
annulus (2.63" - 7.45"). Next, we identify all stars with 
F814W magnitudes brighter than M\\ m = — 3 within the 
two measurement regions. This removes the flux from 
the six brightest stars in Figure [8] We subtract the flux 
contribution from these bright stars from the cluster and 
background flux totals. Finally, we subtract the remain- 
ing background flux (representing the summed light of 
field stars with F814W >My aa ) and obtain our final un- 
resolved fluxes for the cluster. 



The cluster's properties were determined through % 2 
minimization using the unresolved flux values. For 
PC1017, the unresolved light analysis gives a best fit age 
of 40 Myr, with log(M c/ ) = 3.62 M and A v = 0.42 mag. 
All models with x 2 < 1 ± (1 / y/n) were considered good 
fits, where n represents the degrees of freedom, 3 in this 
case. From the distribution of good fits, the 68% confi- 
dence interval gives an age range of 16 - 64 Myr, a mass 
range of log(M c /) = 3.56 - 3.67 M Q , and an extinction 
between 0.20 and 0.64 mag. These values are consistent 
with both the spectroscopically determined properties, 
as well as those predicted by the discrete models. 

We have also analyzed this cluster using traditional in- 
tegrated light fitting. This method gives a best fit age 
of 80 Myr, log(M c ,) of 3.86 M , and A v of 0.28 mag. 
The 68% confidence intervals were 64 - 96 Myr, with a 
log(M c /) between 3.81 and 3.90 Mq, and an extinction 
between 0.17 and 0.39 mag. Fitting the total integrated 
light therefore gives a derived age that is too old, due 
to the presence of the luminous red giant. In contrast, 
the unresolved flux results in an age that is more consis- 
tent with the ages determined spectroscopically, by the 
discrete models, and from isochrone fitting. 

We find that the mass estimated by the unresolved 
method is lower by ~ 0.2 dex compared to the other 
three methods, which is of the same order as the bias 
towards low masses found in S|4l However, the unresolved 
mass estimate is still within the error bars determined 
by the isochrone fitting as well. The underestimate of 
the cluster's mass could indicate that the mass-to-light 
scaling needs to be recalibrated for use with unresolved 
light. 

Though there are differences in the derived ages, 
masses, and extinctions for the different methods, the 
properties determined by the unresolved method agree 
well with those determined by the independent estimates. 
The differences in extinction can mostly be attributed 
to the age-extinction degeneracy. The age determined 
by isochrone fitting was the youngest, and it also had 
the highest extinction. Alternately, this cluster could be 
slightly older, with a lower extinction, as estimated by 
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the discrete and unresolved models. 

6. DISCUSSION 

6.1. Applicability 

Analyzing the unresolved component of a cluster's light 
shows promise for reducing some of the stochastic effects 
associated with deriving the properties of star clusters. 
This method is best suited for low mass clusters, where 
stochastic issues make normal integrated light methods 
fail, and where there are small numbers of bright stars 
that we wish to subtract. The dispersion is still signifi- 
cant for clusters of 10 3 M , but this method can produce 
more accurate results for these lower mass clusters than 
traditional integrated light fitting. It can be used for 
clusters in galaxies from ~ 1-3 Mpc, where full CMD fit- 
ting is not possible, but a small number of bright stars 
can be individually resolved. The value of Mn m should 
be optimized for the targeted clusters' age range. 

6.2. Best Choice of Mu m 

The choice of cutoff magnitude M\ im was based on try- 
ing to achieve stability in the unresolved component of 
the clusters' light, while preserving a strong correlation 
of flux with age. The values of M\i m discussed in this 
paper were chosen to work well with the PHAT data set. 
Several things should be kept in mind when trying to 
decide what value of M lim should be used. At distances 
greater than a few Mpc, even HST imaging will make 
it difficult to resolve stars down to moderately bright 
magnitudes (e.g., M 8 i4 — —3). However, in these cases, 
the method should be adjusted to use a brighter Mn m at 
the expense of limiting the age range one can probe. A 
brighter cutoff would be better suited to younger and/or 
less extincted clusters, while a fainter cutoff would yield 
better results for older and/or more extincted clusters. 

One difficulty associated with subtracting the resolved 
star fluxes is that completeness of resolved photometry 
can vary greatly as a function of radius within the cluster. 
Miim should be above the completeness limit of the data 
since accurate stellar photometry of the resolved stars 
is needed to subtract off their flux. If My lm is bright 
enough and the number of bright stars we wish to sub- 
tract is small enough, then resolving these few bright 
objects should be possible in most cases. This will be 
the case for lower mass clusters, for which this method 
is most useful. To optimize this method for a general 
use, Mii m could be a parameter that is solved for during 
the x 2 minimization as well. This would allow the value 
of Mi; m to be optimized for the age of each cluster, and 
to reduce stochastic effects for older clusters as well as 
younger clusters. 

6.3. Field/Foreground Contamination 

Another benefit to this method not previously dis- 
cussed in that it minimizes bright field/foreground star 
contamination. As long as the star in question is brighter 
than the chosen Mi; m , its flux is not included in the anal- 
ysis. Therefore, the determination of cluster membership 
does not affect the derived cluster's properties. 

6.4. Model Uncertainties 

Our analysis depends on the accuracy of the models 
used and the assumption that dust attenuation acts in 



a predictable way. These factors may cause the actual 
uncertainties in derived properties to be higher. One ad- 
vantage of looking at the unresolved light is that this will 
be less sensitive to model uncertainties for massive post- 
main sequence stars, which are quite substantial. One 
potential complication, however, is that the separation 
between main sequence and post-main seq uence stars is 
not a lways as clean as the models predict (Larsen et al. 
2011 1, which again complicates the choice of Mn m . 



6.5. Future Work 

Eventually the best way to study a cluster is to com- 
bine the study of the unresolved and resolved portions 
of its light. The resolved stars can be analyzed using 
isochrone fitting, and the unresolved color could be ana- 
lyzed separately. Then the results from these two meth- 
ods can be compared to ensure consistency. 

A more complete study of a variety of clusters is needed 
to further show the applicability of this method. This 
would include optimizing the value of M\{ m to be solved 
for during the fitting process, investigating the effects of 
crowding and blending, further comparisons with other 
age-dating methods, and extending this method out to 
larger distances. 

7. CONCLUSIONS 

We investigated the properties of partially resolved 
stellar clusters using simulations and found that stochas- 
tic variations in color can be greatly decreased when con- 
sidering only flux below a limiting magnitude Mu m . This 
unresolved light component utilizes the stability of the 
main sequence light along with the age information in the 
main sequence turnoff but eliminates the often stochasti- 
cally sampled upper end of the stellar mass function. By 
using only this unresolved component of the flux, we have 
shown that we can derive accurate age, mass, and extinc- 
tion determinations with a variety of simulated clusters. 
The improvements over traditional integrated light fit- 
ting are most evident for lower mass clusters where the 
effects of stochasticity are greatest. This method was 
also applied to an M31 cluster, and results using the 
unresolved method were comparable to the properties 
determined spectroscopically, from discrete models, and 
isochrone fitting. 

This new technique will be especially useful for lower 
mass clusters (less than a few 10 4 M ), crowded clusters, 
and clusters in nearby galaxies with only a few resolved 
stars. This method can also be used as a sanity check 
for clusters whose age and mass determinations come 
from isochrone fitting. It also allows for the potential 
of obtaining reliable property determination for clusters 
too far away to be analyzed with isochrone fitting meth- 
ods. In this situation, the unresolved technique provides 
greater accuracy while still utilizing continuous models, 
without needing spectroscopy or computationally inten- 
sive discrete models. 
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